Finite volume effects in Bk with improved staggered fermions 
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We extend our recent unquenched (Nf = 2+1 flavor) calculation of Bk using improved staggered 
fermions by including in the fits the finite volume shift predicted by one-loop staggered chiral 
perturbation theory. The net result is to lower the result in the continuum limit by 0.6%. This shift 
is slightly smaller than our previous estimate of finite volume effects based on a direct comparison 
between different volumes. 
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We have recently reported a result for Bk using im- 
proved staggered fermions with all errors controlled and 
with a total error of 6% (We refer to this paper as 
SWME in the following; a recent update including an 
additional lattice spacing is given in Ref. 0.) Although 
the dominant error in this result is from uncertainty in 
the matching factors, an important subdominant source 
of error is that arising from our use of a finite volume. In 
SWME we estimated this error to be 0.85% by comparing 
the result obtained on two lattices of different volumes 
(ensembles C3 and C3-2, as discussed below). Here we 
revisit the finite volume (FV) error using an alternative 
approach: repeating our chiral fits using forms predicted 
by one-loop SU(2) staggered chiral perturbation theory 
(SChPT) but now including FV effects from pion loops. 

We did not previously include one-loop FV effects be- 
cause we found that fitting with the FV form using CPUs 
on our extensive dataset took too long (of order two 
months). We now use GPUs (Graphics Processing Units) 
which allows us to reduce the fitting time down to a 
few hours, while maintaining double-precision through- 
out. Thus we can undertake the many analyses needed 
to estimate systematic errors. 

To extrapolate our data to the physical light quark 
masses and to the continuum limit we use the functional 
form predicted by SChPT d ||. Specifically, our cen- 
tral result in SWME used partial next-to-next-to-leading 
order (NNLO) SU(2) SChPT fits, and we consider here 
only such fits. The key extrapolation is in the valence 
down-quark mass, m Xl or more precisely in the mass- 
squared of the valence pion with flavor xx, which we call 
Xp. The NNLO form contains 3 parameters (at a fixed 
lattice spacing a and fixed sea-quark masses): the coeffi- 
cients of (i) the constant + chiral logarithm term, (ii) the 
term linear in Xp and (iii) the term quadratic in Xp. For 
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details see Eqs. (41-45) of SWME. Finite volume correc- 
tions only impact the chiral logarithms. In SWME these 
are written in terms of the functions: 



e(x) = x[io g (x/^) + s F 1 v (x)} , 
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where /xdr is the scale introduced by dimensional regu- 
larization. The FV parts, left out in SWME but included 
here, are 
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where L is the spatial box size and n = (ni, ria, Tia, n^) is 
a vector of integers labeling image positions. The norm 
In I is defined as 
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with Lt is the (Euclidean) temporal box size. Kq and K\ 
are the standard modified Bessel functions of the second 
kind, which fall exponentially for large x. 

The expressions (|3|) and ((4} for the FV corrections 
cease to be valid when ML < 1. One then moves from the 
so-called "p-regime" into the "e-regime" , and a different 
power-counting applies. Our calculations are done with 
valence pions satisfying ML > 3, which, as the following 
results indicate, appears to be above the minimum value 
at which it is appropriate to use the above expressions. 

The lattice ensembles used in this paper are those used 
in SWME to do the continuum extrapolation, namely 
(in the notation of SWME) C3 (coarse), Fl (fine), SI 
(superfine), together with a new ensemble Ul (ultrafine). 
The latter has a nominal lattice spacing of 0.045 fm, sea 
quark masses am? — 0.0028 and am s = 0.014, and size 
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64 3 x 192. Our results for this fourth lattice spacing are 
preliminary, based on only 305 configurations. We also 
use the larger volume coarse ensemble, C3-2. 

On each ensemble we calculate Bk using 10 differ- 
ent sea quark masses, ranging down from rs mf ys to 

wP hys /10. The precise values are given in SWME for 
the SI, Fl, C3 and C3-2 ensembles, while for the Ul en- 
semble they are 0.0014 x n with n = 1, 2, 3, . . . , 10. For 
the chiral extrapolation we use the lightest four valence 
masses for m x , while for the valence strange quark (which 
we call m y ) we use the three heaviest valence masses to 
extrapolate to the physical strange-quark mass. We do 
the fit in two stages, first (the "X fit") extrapolating in 
the valence d quark mass to the physical value while hold- 
ing the valence s quark mass fixed, and second (the "Y 
fit" ) extrapolating linearly to the physical valence s quark 
mass. Altogether this procedure is labeled a "4X3Y" fit. 

To calculate the finite volume corrections Sj v and 5f v , 
we must truncate the sums in Eqs. §3§ and (U). We keep 
image vectors satisfying \n\ < n max , with n max deter- 
mined separately for each sum and for each value of ML, 
using the following criterion. For \n\ < n max , a shell in 
image space of radius |n| and unit thickness should give 
a contribution larger than the desired precision times the 
leading (|n| = 1) term. Thus we keep all n satisfying 

Ki(\n\ML) 



[4tt 



> e x [6Kx (ML)] for Sf v , (6) 



[4ttH 2 ] x K Q (\n\ML) > e x [QK (ML)} for (5f v , (7) 

with e = 1.0 x 10~ 14 since we want double-precision ac- 
curacy. Note that on the left-hand side we are assuming 
that the dominant contribution is from spatial images, 
which is a reasonable approximation. 

The values of n max depend, of course, on ML. For v 
we find, for example, that for ML = 0.4, 2.0 and 4.0 that 
"■max = 93, 18 and 10, respectively. The values for 8f v 
are slightly larger, although comparable. 

Using a single core of the Intel i7 920 CPU, it takes 
about two months to run the analysis code on our full 
dataset in double precision including FV corrections, at 
a sustained speed of about 0.5 Gigaflop/s. By contrast, 
using the Nvidia GTX 480 GPU, and optimizing our 
code using the CUDA library, we obtain a sustained per- 
formance in double precision of 64.3 Gigaflop/s, 38% of 
peak. 

We first compare the results of fitting with and with- 
out FV corrections on the C3 ensemble. We note that 
on this ensemble (and on C3-2) we have about 9 times 
the number of measurements as on the finer ensembles 
and thus considerably smaller statistical errors. Figure [Tj 
shows the X-fits for our heaviest valence strange quark 
(am™ 1 = 0.05). Although our lightest valence pion mass 
corresponds to ML = 2.7 on these lattices, the difference 
between the infinite and finite volume fits is not signifi- 
cant until much smaller values of ML. This is in apparent 
contradiction with rule-of-thumb that one has large FV 
effects when ML < 3. This conundrum is resolved by the 
presence of taste breaking. Almost all the pions which 
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FIG. 1. Bk(1/cl) vs. Xp on the C3 ensemble for am y = 0.05. 
The fit type is 4X3Y-NNLO. The dashed [red] line shows the 
result of fitting to the infinite volume form, while the solid 
[blue] line shows the fit with FV corrections included. The 
[red] diamond and [blue] square correspond to the "physical" 
Bk values obtained (as described in the text) from the infinite 
and finite volume fits respectively. 



appear in the loops have non-Goldstone tastes, and thus 
are heavier than the Goldstone pion by shifts of 0(a 2 ). 
These shifts push ML up to values considerably larger 
than 3, except for a small contribution from the Gold- 
stone pion itself. 

The net result is that fitting with the FV form leads to 
a very small shift in the extrapolated "physical" values 
of Bk , as shown in the first row of Table HI These values 
are obtained by using the fit function to (a) extrapolate 
to the physical valence d mass, (b) set all taste-splittings 
to zero, (c) set the light sea-quark pion mass-squared to 
its physical value, and (d), in the case of the FV fit, set 
the volume to infinity. We use quotes around "physical" 
because these values of Bk have yet to be extrapolated 
to the physical strange quark mass and to the continuum 
limit. Note also that Bk is here matched to a continuum 
operator renormalized at a scale 1/a, so that results from 
different lattice spacings are not directly comparable. 



TABLE I. "Physical" Bk (NDR, 1/a) obtained as described in 
the text from fits without and with finite volume corrections, 
and the percentage difference between the two. The valence 
strange quark mass is set to to its heaviest value on each 
ensemble. Errors are statistical. 
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B K (FY) 


AB K 


C3 


0.5734(46) 


0.5743(46) 


+0.159(2)% 


C3-2 


0.5784(46) 


0.5785(46) 


+0.0319(3)% 


Fl 


0.5225(111) 


0.5199(110) 


-0.505(14)% 


SI 


0.4914(65) 


0.4898(65) 


-0.329(6)% 


Ul 


0.4780(92) 


0.4757(92) 


-0.474(14)% 
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From the Table, we see that the FV fit leads to a result 
that is 0.16% higher on the C3 ensemble. Although this 
difference is much smaller than the statistical errors in 
the individual results, it is statistically highly significant. 
This is possible because of the high degree of correlation 
between the infinite volume and FV fits. The same holds 
true on the other ensembles. 

Figure [5] shows the corresponding fits on ensemble of 
larger lattices, C3-2. As expected, the difference between 
the infinite and finite volume curves does not become sig- 
nificant until smaller values of the pion masses than in 
Fig. [TJ The "physical" Bk values resulting from these 
two fits are essentially the same. This supports our claim 
in SWME-1 that we can essentially treat the C3-2 ensem- 
ble as having infinite volume when calculating Bk- 
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FIG. 2. Bk(1/o) on the C3-2 ensemble. Notation as in 

Fig.m 
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FIG. 3. Sjr(l/a) on the Fl ensemble. Notation as in Fig. [T] 
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A natural question is whether the FV fit leads to bet- 
ter agreement between the C3 and C3-2 ensembles. The 
answer is a qualified yes. The FV shift does bring the 
two results closer, although it removes only 20% of the 
difference. It should be noted, however, that the results 
on these two ensembles, which are statistically indepen- 
dent, agree within 1-cr already before the FV correction. 
Thus the improvement is marginal. 

Figures [21 HI and [5] show the fits with and without FV 
corrections on the Fl, SI and Ul ensembles, respectively. 
We see that the FV corrections are of opposite sign to 
those on the C3 and C3-2 ensembles, that FV effects set 
in at a larger value of Xp as one reduces a, and that FV 
effects at the physical value of Xp grow rapidly as a is 
reduced. All these changes are due to the reduction in the 
size of taste breaking as a is decreased, so that the values 
of ML for non-Goldstone pions are much closer to those 
of the Goldstone pion, enhancing FV effects. It turns out 
that the enhancement is greater for terms contributing 
positively than those contributing negatively, leading to 
the change in sign. It is clear from the curves that our 



FIG. 4. Bx{i/a) on the SI ensemble. Notation as in Fig. [T] 
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FIG. 5. B K (l/a) on the Ul ensemble. Notation as in Fig.[T] 
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lightest valence pions are close to the smallest values that 
can be used without encountering large FV corrections. 

The corresponding values of "physical" Bk are given 
in Table|H The FV shifts are now larger, as large as 0.5% 
in magnitude, although still small compared to other sys- 
tematic errors. We stress that one cannot gauge the size 
of the FV shift at the physical value of Xp by looking 
at the difference between the two curves for this value. 
This is because the fits are being done at larger values of 
Xp, where the FV corrections are smaller. 

In Fig. [SI we show how the continuum limit is impacted 
by the inclusion of FV corrections. Here we plot Bk after 
extrapolation to the physical valence strange quark mass 
(using 3Y fits), and after running in the continuum to a 
common scale of 2 GeV. The overall effect of fitting with 
FV corrections is to reduce the continuum result from 
0.5260(73) to 0.5229(72), i.e. by 0.59 ± 0.01%. 
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inclusion of FV effects arising from pion loop diagrams. 
The net result is that, after continuum extrapolation, we 
find that including the FV shifts leads to an (infinite vol- 
ume) result for Bk which is 0.59% smaller than that ob- 
tained when fitting with the infinite volume forms. The 
smallness of this shift confirms that our values of ML 
(which range down to f=a 2.7) are (barely) large enough. 
In this regard, we are helped by the presence of taste 
breaking, so that most of the pions appearing in the loops 
are heavier than our lightest pion. 

As is well known, one-loop ChPT forms can under- 
estimate the size of FV shifts, due to the contributions 
from higher-order terms. This underestimate can be by 
as much as a factor of 2 for our range of ML (sec, e.g., 
Ref. [5]). Thus, our new procedure for accounting for FV 
errors is to take the central value from the fits includ- 
ing FV effects, and take the systematic error to be the 
size of the FV shift, namely 0.6%. This estimate is a 
little smaller than our previous estimate of the FV error, 
0.89%, which was based on the difference between the re- 
sults on the C3 and C3-2 ensembles. We think that our 
new procedure is more reliable, since it includes a con- 
tinuum extrapolation, and is based on the theoretically 
expected functional form (which, in particular, includes 
the expected dependence on a). Our previous estimate 
also did not account for the possibility that the difference 
on the two ensembles C3 and C3-2 could be statistical. 
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FIG. 6. Continuum extrapolation of Sjr(NDR, 2 GeV). [Red] 
octagons show results obtained using SU(2) 4X3Y-NNLO 
SChPT fits without the finite volume corrections, while [blue] 
crosses show the results obtained with FV fits. The points at 
a = are obtained by linear extrapolation in a 2 . 

In conclusion, using CPUs we have been able to im- 
prove the chiral extrapolation of our data for Bk by the 



C. Jung is supported by the US DOE under con- 
tract DE-AC02-98CH10886. The research of W. Lee 
is supported by the Creative Research Initiatives pro- 
gram (3348-20090015) of the NRF grant funded by the 
Korean government (MEST). The work of S. Sharpe is 
supported in part by the US DOE grant no. DE-FG02- 
96ER40956. Computations for this work were carried out 
in part on QCDOC computers of the USQCD Collabora- 
tion at Brookhaven National Laboratory. The USQCD 
Collaboration are funded by the Office of Science of the 
U.S. Department of Energy. 



[1] T. Bae, Y.-C. Jang, C. Jung, H.-J. Kim, J. Kim, J. Kim, 

K. Kim, W. Lee, S. R. Sharpe, and B. Yoon, Phys.Rev. 

D82, 114509 (2010), arXiv: 1008.5179. 
[2] T. B. Boram Yoon, Y.-C. Jang, H.-J. Kim, J. Kim, J. Kim, 

K. Kim, W. Lee, C. Jung, and S. R. Sharpe, PoS (Lattice 

2010), 319 (2010), arXiv:1010.4778. 



[3] W. Lee and S. R. Sharpe, Phys. Rev. D60, 114503 (1999), 

hep-lat/9905023. 
[4] C. Aubin and C. Bernard, Phys. Rev. D 68, 034014 (2003), 

arXiv:hep-lat/0304014. 
[5] G. Colangelo, S. Durr, and C. Haefeli, Nucl.Phys. B721, 

136 (2005), hep-lat/0503014. 



